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Abstract 



This paper presents a comprehensive study of the functions Vmix) = 
P(^_^_l) J^it^ — xP)'^e~^^ dt for X > 0, m > —1 and p > 0. For large x these 
functions approximate x^~^. The case p = 2 is of particular importance 
because the functions V^{x) ~ 1/x can be regarded as one-dimensional 
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regularizations of the Coulomb potential l/|x| which are finite at the origin 
for m > — 2- 

The limiting behavior and monotonicity properties of these functions 
are discussed in terms of their dependence on m and p as well as x. Several 
classes of inequalities, some of which provide tight bounds, are established. 
Some differential equations and recursion relations satisfied by these func- 
tions are given. The recursion relations give rise to two classes of polynomi- 
als, one of which is related to confluent hypergeometric functions. Finally, 
it is shown that, for integer m, the function 1/V^{x) is convex in x and this 
implies an analogue of the triangle inequality. Some comments are made 
about the range of p and m to which this convexity result can be extended 
and several related questions are raised. 
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1 Introduction 

1.1 Definitions and background 

In this paper we study the functions 



2e^'' 



fOO 

Vm{x) = TT7^^-—T I {f-x're-'"dt, m>-l (1) 



V-iix) = 
and their generahzations, 



r(m + l] 

1 



x\ 



K^(^) = f(£^/ (t'-^Te-^'dt (2) 



a-P 



for < p < oo. These functions are well-defined for x > and can be extended 
to complex m with ^{m) > — 1. For ^{m) > — | they are also well-defined for 
X = 0. Using symmetry or the equivalent forms (^ and (0) below, they can be 
extended to even functions on R or R\ {0}. However, it suffices to consider only 
non- negative x and in this paper we restrict ourselves to that. We also restrict 
ourselves to real m. 

Letting p = 2 in (^ yields (||). However, because this case is more important 
in applications we often drop the superscript and and simply write Vm{x) for 
V^ix). 

Our interest was motivated by studies of atoms in magnetic fields where 
these functions arise naturally for integer m. Vm can be regarded as a (two- 
dimensional) expectation of the (three-dimensional) Coulomb potential l/|r| with 
the state 7m(r, ^) = ^ ^ ^-ime^m^-r /2 (^-^j-^g^-g -^g have used cylindrical coordi- 
nates r = (x, r, 6) with the non-standard convention r = ^Jy'^ + z^ if r = (x, y, z) 
in rectangular coordinates). The state 7^, describes an electron in the lowest of 
the so-called "Landau levels" with angular momentum m in the direction of the 
field. In this context, it is natural to rewrite (|ID in the form 



00 2m -r^ 



V^(x) = — / rdr, (3) 



00 ^m^-u 



r(m+i)7o V^^ 



u 



■du (4) 



for m > —1. In this form, it is easy to see that Vm{x) ~ ^/\x\ for large x. The 



importance of Vm goes back at least to Schiff and Snyder |T£] and played an 
essential role in the Avron, Herbst and Simon study of the energy asymptotics 
of hydrogen in a strong magnetic field. More recently work in astrophysics and 



the work of Lieb, Solovej and Yngvason [|I^ on asymptotics of many-electron 
atoms in strong magnetic fields has renewed interest in this subject. Motivated 
by the LSY work, Brummelhuis and Ruskai 0, H have developed one-dimensional 
models of many-electron atoms in strong magnetic fields using the functions Vm{x) 
as one-dimensional analogues of the Coulomb potential. 

In the case of many-electron atoms, the anti-symmetry required by the Pauli 
exclusion principle suggests replacing the simple "one-electron" expectation above 
by an N-electron analogue in which the state 7^ is replaced by a Slater determi- 
nant of such states. This is discussed in detail in [Q where it is shown that, in 
the simple case corresponding to m = . . . A^ — 1, the analogous one-dimensional 
potentials have the form 



VAx) 



1 ^-^ 



m=Q 



In Section |^ we obtain recursion relations for Vm which, in addition to being of 
considerable interest in their own right, are extremely useful for studying poten- 
tials of the form (^. 

For m = the function 



occurs in many other contexts and is sometimes called the "Mills ratio" pA 
Although it has been extensively studied, the class of inequalities we consider in 
Section ^?1| appears to be new (although some of our bounds coincide with known 
inequalities in other classes) and the realization that l/Vo(x) is convex seems to 
be relatively recent []T^ |T^, 



The replacement of x by x^ in (|^) has been considered by Gautschi [jTO| and 



Mascioni |]T2|, who (after seeing the preprint |T^) extended the results of Section 



4?1| to this situation. 

For the analysis of this generalization, it is useful to observe that (|^) can be 
rewritten as 






In this form, it is easy to verify that VJ^{x) ~ x^ ^ for large x and that for p = 1, 
V},(x) = 1 for all m . 



Our first result shows that V^{x) is continuous in m and that our definition 
for m = — 1 is natural. 



Proposition! For all x > 0, lim z^ ^V^i 



X] 



-1+ 



Proof: Note that 



X- 



p-i 



p-i 
(x'P + u) p 



i^+^y 



— so that (^ implies 



l-xP-^VP{ 



x] 



T{m+l)Jo 



M™e-" 






du. 



Since r(m + 1) becomes infinite as m ^ —1, the desired result follows if the 
integral on the right above remains finite. To see that this is true, it is convenient 

to let g{z) = - (l — {I + z)~p~ ] and note that (||) implies 



X- 



p-i 



VZ(x) < 



1 



T(m + l) 



oo ^m+lg-u 



XP 



9 



u 



X'P 



du 



(9) 



It is easy to see that the large u portion of this integral causes no problems 
since \g{z)\ is bounded by a polynomial in z when z > 1. (For p > 1, it is 
bounded uniformly by 1; for < p < 1, it is bounded by a polynomial, namely 

\9{.z)\ <! + (! + zY , where /;; G N , A; > -.) To see that it is also well-behaved 

for small m, we first note that for p > 0, (1 + 2;) p is analytic for '^{z) > —1. 
Then g{z) has a removable singularity at 2; = and can be extended to an analytic 
function on "^{z) > — 1. Thus, for small m, the integrand behaves like a;~^M'""'"^e~'" 
which ensures that the integral in (^ is finite for m = —1. QED 

The rest of this paper is organized as follows. In the next part of this section we 
summarize the properties of Vm in the important case p = 2. We then conclude the 
Introduction with a summary of convexity results, including some open questions. 
In Section | we state and prove the basic properties of V^ for general p. In Section 
^ we derive recursion relations for V^ and study their consequences. Among these 



is a connection with confluent hypergeometric functions. In Section |4.1| we prove 
some optimal bounds for Vq. The optimal upper bound had been established 
earlier independently by Wirth [|l^ and by Szarek and Werner [jl^ who also 



showed that the upper bound is equivalent to the convexity of I/Vq. In Section 



]2| we discuss several classes of inequalities, beginning with optimal bounds on 
Vo{x). We then consider optimal bounds on the ratio Rm{x) = Vm{x)/Vm-i{x) 
and show that these have important consequences. In particular, we show that 
the upper bound is equivalent to the convexity (in x) of 1/Vm{x) and that the 
ratios increase with x. Proofs of the ratio bounds are then given in Section ^ 
where we also consider extensions to other p. Because the proof of the ratio 
bounds is via induction on m, the results of Sections 6 and 7 are only established 
for integer m. However, we believe that they hold for all m > — 1. 

1.2 Properties of Vm{x) 

We now summarize some properties of Vm{x) along with comments about the his- 
tory and brief remarks about the proofs. Unless otherwise stated, these properties 
hold for m > —1 and x > 0. 

1 1 

> Vrr,ix) > 



\/x^ + m \/x^ + rra + 1 

where the first inequality holds for ?7i > and the second for m > — 1. 

To prove the upper bound, which appears to be new, observe that yU = 
[M'"~^e~"/r(m)](iti is a probability measure on (0, cxd). For fixed x, one can 
then apply Jensen's inequality to the concave function /^^(m) = u{u+x^)~^l'^ 
to obtain 



mJo m \Jq 



r(m) 



—fx{m) = -^== 
m Vx^ + 



m 



The lower bound was proved earlier (at least for integer m) by Avron, Herbst 
and Simon p[ who applied a similar argument to the probability measure 

[M'"e~"/r(r7i + l)]^^ and the convex function fx{u) = (-u + x^)^^/^. 

b) Vm{x) is decreasing in m. In particular, Kn+i(x) < Kn(x) < — . 

X 

The first inequality follows easily from property (a) which implies 
Vm{,x) < -j= — < Vm-i{x). Alternatively, one could use integration by 
parts on @). The second inequality is easily verified from the integral rep- 
resentation (H). That Vm{x) also decreases with m for non-integer jumps 
is more difficult, and the proof is postponed to Section || where it follows 
from the more general Theorem ^. 

6 



c) The expression mVmix) is increasing in m > — 1, m G R. 

For integer jumps this holds for m > —1. Indeed, it is obvious that — V^-i < 
■ V^o < Vi- For integer jumps with m > 1, one can use property (a) to see 
that 



mVr„(x) > 



m 



> 



m — 1 



y/x^ + 771+1 y/x"^ + 777—1 



>(m-l)Kr,_i(x). (10) 



The proof for general 777 is postponed to Theorem ^ in Section |^. The 
fact that Vm{x) is decreasing in 777, while mVmix) is increasing gives an 
indication of the delicate behavior of Vm- 

d) For 77i > —1/2, the definition of Vm{x) can be extended to x = and 



K^(O) 



For integer 777, this becomes 

(2m)! 



Vm{0) 



n 



r(777 + |) 
r(777 + l)' 

l-3-5...(2m- 1) 



'11^ 



22m(^^!)2 V 2 -4 -6... (2m) 

while for large m Stirling's formula implies 

1 



(12) 



K^(0) 



77i 



I \ m 
2 



771 



e\i/2 

777/ 



777 



(13) 



which is consistent with property (a). Boyd |^, T^\ has proved the more 
precise estimates 



'^ + 4 ~'~ 32m+48 



m + 



< vuo) < 



1 



"^ "•" 4 "•" 32m+32 



(14) 



e) For all m > 0, Vm satisfies the differential equation 

V^^ix) = 2x {V„, - Vm-l) ■ 
This can easily be verified using integration by parts in (^). 

f) For each fixed 777 > 0, Vm{x) is decreasing in x. 
This follows directly from (b) and (e). 



(15) 



g) For a > 0, the expression aVm{ax) increases with a. Hence aVm{ax) > V{x) 
when a > 1 and aVm{ax) < V{x) when a < 1. 

This property follows easily from the definition (^) or (Q) and the observa- 
tion that = = — , is increasing in a. It is used in the proof 
Va2x2 + u v^S^T^ 

of Theorem ^ and is important in the study of one-dimensional models for 
atoms in magnetic fields in which the electron-electron interaction takes the 
form of convex combinations of 4=Kn f \^j ^k\ \ ^ 

h) Vo{x) is convex in a; > 0; however, Vm{x) is not convex when m > |. 

For m = 0, the differential equation (|15D becomes Vq{x) = 2[xVo — 1]. Since 
xVq = r„°° / ^ du is increasing for x > 0, it follows that Vo(x) is convex. 

When m > I it follows from (|^) and (b) that lim^^o Kn(^) — 0- Since V^ 
is negative, V^ must decrease on some small interval (0,a;o). One can also 
show lima-^oo Kra(^) — 0' so that one expects that there is an xi such that 
Vm is concave on (0, xi) and convex on {xi, cx)). In Section ^ we will see that 
the convexity is recovered for the averaged potential V^ 



rav 
m ■ 



i) For integer m, 1/Vm{x) is convex in x > 0. 

This will be proved in Section [4.2| as Theorem ^ For large x, 1/Kn(x) ~ x 



so that the deviation from linearity is very small and the second derivative 
close to zero. This makes the proof quite delicate and lengthy. 

The convexity of 1/Vm{x) can be rewritten as 

1 1 1 



Using property (g) with a, = |, one easily finds that the convexity of 1/Kn(x) 
implies 



1 1 1 



This subadditivity inequality plays the role of the triangle inequality in 
applications. (See, e.g. 0].) 

j) Asymptotic estimates: 

For large x, it follows from property (a) that 

m 1 ^^ , , m + 1 . , 

< - - V^x) < — — (16) 



2(a;2 + m)3/2 - x ^ ' 2x^ 
8 



The asymptotic expansion 

can easily be obtained from (^). For details let p = 2 in the proof of 
Proposition |^ which gives a similar expansion for p > 1. 

It follows from properties (a) and (c) that Kn(x) decreases monotonically to 
zero for each fixed x as m —>■ oo. In fact, since Vm{x) is decreasing in x for 
all m, it suffices to show this for x = which is easy since Kn(0) < m~^/^. 

k) The Fourier transform is given by 

^ 1 fXi Am+l poo m p-s 

V^{i) = ^ Kn(x)e-'"«(ix = ^= / -—. T^-rds (18) 

^^' V2^J-oo V2^Jo (|eP + 4s)™+i ^ ' 



This follows from (|]) and the standard formula (e.g., see (v) on p. 131 of 
ill) ^1 , ^ =)(^) = A^ /""le-iH"l'l«l'/^)t/s after achange 
in the order of integration. 

1.3 Convexity Summary 

For large x, V^{x) ^ x^~p which is convex in x for p > 1 and concave for 
p < 1. For m > — I these convexity properties can not be extended to V^^^x) 
on all of (0, oo); they would be inconsistent with the differential equation and 
monotonicity properties in Proposition ^ However, as discussed after Proposition 
[TT], the averaged potentials V^^^ have the same convexity as x^~^ on the half-line. 
The convexity of 1/V^{x) is the motivation for Sections ^ and |^. This question 
is already delicate for p = 2 and its verification becomes increasingly difficult for 



larger p. Although, as discussed in Section p^ , we have evidence that convexity 
holds for all p > 2, our methods give this result only for a limited range of p. 
Moreover, because our proof is inductive, we have established convexity and ratio 
bounds of Section [4.2| only for integer m. It would be interesting to find another 
approach which would extend these results to non-integer m and all p > 2. Since 
1/V^{x) ~ x^^^ which is concave for 1 < p < 2, we can not expect convexity of 
1/V^ in this range. 

As discussed above, one important consequence of the convexity of 1/Vm{x) is 
an analogue of the triangle inequality. For all values of p we have [^^(s;)]^^ ~ x 



for large x which suggests a triangle inequality of the form 

It would be interesting to know the range of p (and m) for which this holds. For 
p > 2 the convexity of 1/V^{x) implies only the weaker inequality 

Finally, one could also ask if Kn(x) is convex in m. In particular, is 2V^(x) < 
Vm+i{x) + Vm-i{x) or, equivalently by (|T5|) , is V^{x) increasing in m? 

2 General p 

We now study the basic properties of V^ in detail. As one would expect from 
V^(x) ~ x^~^, the behavior of V^ is often quite different for p > 1 and p < 
1. At the boundary, p = 1, V^{x) = 1 for all x. Proposition |^ describes the 
monotonicity and limiting behavior of V^{x) as p varies with m and x fixed. 
Proposition || gives a simple expression for V^ in the special case that 1/p is an 
integer. 

The next four results generalize properties of Section [L^ to general p. Propo- 
sition ^ generalizes the inequalities from property (a); Proposition ^ generalizes 
properties (d), (e), (f), and (g); and. Theorem |] extends the monotonicity prop- 
erties (b) and (c). Moreover, the proof of monotonicity for non-integer jumps is 
provided here. Finally, Proposition |^ gives the asymptotic behavior of V^{x) for 
large x when p > 1. 

Proposition 2 Let m > —1 and x > be fixed. Then 
(i) limp_oV;^(a;) = oo. 
(ii) For all X > 1, V^ is decreasing in p. Moreover, 

if X > 1, lim V^{x) = 0, and 

p— >oo 

^fx = l, limKi:(l) = fT^/ ^V^- 
p^oo l{m + lj Jq 1 + u 

(iii) For all < x < 1 and m > 0, limp^oo Kn(^) = m- 



10 



Proof: We use the expression (0) for V^. 
(i) Since linip^ol^;^ + m)^/^ = oo for x > 1, 

lim Vl{x) > lim — / ^ ^ = oo. 

(ii) Differentiating (|^) yields 

dp r(m + i)yo p2(xp + M)^ 



(1 - p)x^ In(xP) - (x^ + u) ln(x^ + m) 



du. 



For X = 1 or p = 1, the first term in square brackets above is zero leaving a 
quantity which is clearly negative. When both x > 1 and p > 1 both terms in 
([19| ) are clearly negative. When x > 1 and < p < 1, the quantity in square 
brackets in (p!9D is negative since 

(1 - p)xP In(xP) - (xP + u) ln(xP + u) 

< xP In(xP) - (x^ + u) \n{xP + u)<0. 

The last inequality follows from the fact that the function f{w) = wlnw is 
increasing for w > 1. Thus, j-Kn(a;) < for all x > 1 and for all p G (0, oo). 

(iii) Since limp^oo(a;^ + uY^^ = 1 for x < 1, 

m—l„—U' 
, '/: I ^ I 7 ^ / 



hm V^{x) = / n-^e-"rfw = -.QED 

p^oo i (m + 1) Jq m 

Remark: The behavior for x < 1 depends upon m. For "small" m, there is a po 
such that V^^ decreases (below 1) on (0,po) and increases on (po, oo) to ^. For 
"big" m , V^ simply decreases to — as p increases. 

The next result shows that in the special case that 1/p is an integer V^ reduces 
to a polynomial in x^ = x^'" of degree n — 1. 

Proposition 3 For n G N, n > 2, 

for all X > and m > —1. 

11 



Proof: It follows from (|^) that for p = - 

T{m+l)Jo 

When n is an integer > 2, the result then follows easily from the binomial expan- 
sion applied to {x^ + u)"'~^ and the definition of the F-function. 

Proposition 4 For all x > 

^ -<V^{x)< ^—j^ forp>l, 



JJ— 1. IIL \ / J--— ^ 

[x^ + m + 1) p [xP + m) p 

where the first inequality holds for m > —1 and the second for m > 0. 

l~p 1 — p J_ 

{x^ + m + 1) p > V^{x) > {x^ + m) p for - < J9 < 1, 
where the first inequality holds for m > —1 and the second for m > 0. 

V^{x) > {x'' + m + 1)^ for <p< -, 
and the inequality holds for m> —\. 
Proof: The proofs are done using Jensen's inequality as in property (a) of Section 



1.2 



Proposition 5 For all x > 

(i) Form> -i, 1/^(0) zs defined and \/^(0) - ^ 



r(m + l) 
(ii) For all m >0, x > 0, V^ satisfies the differential equation 

^Vlix) = px^-' {V^{x) - V^_,{x)) . (20) 

(iii) For all m > —1, V^ is decreasing in x, if p > I, identically equal to 1 for 
all X, if p = I, increasing in x, if p < 1. 

(iv) Let m > —1 and x > 0. For a > 0, the expression a^^^V^lax) increases in 
a, if p > 1 and decreases in a if p < 1. 

12 



Proof: The proofs are straightforward extensions of those given in Section |1.2| . 
In (iii), one can verify that V^ is also increasing for < p < | by computing the 
derivative directly. 

Theorem 6 For each fixed x > 0, and for m in the region m > —1, 

(i) V^{x) is strictly decreasing in m forp > 1 and strictly increasing in m for 
p <1. 

(ii) mV^{x) is strictly increasing in m for p > 1 and strictly decreasing in m 
for p <1. 



Proof: To prove (i) we differentiate (|^ to get 



—V'Pix) = ^ 

dm "" r(m+l) Jo ixP + uV'v 






(21) 



Using the same procedure as that used (see, e.g. |I|, |n|]) to obtain the standard 
integral representation 



ij{z) 



r(£) 

r(^) 



s s(l + s) 



ds 



(22) 



one finds 
1 



r(m + 1) 7o {x'P + 

1 



— -du 



r(m + ll 



U) p 
°° ds '"^ 



=0 -^ Jo 



\e — e 



h t:x^« 



OOg-S 



Vlix) / ds 



xP + u) p 

ds f°° e-'"w'^dw 



OOg-s 



s T{m+l)Jo s{s + ir^pJo [xP{s + l)+w]^~p 

Vl{x) I —ds - r Vl{x{s + 1)^) ^ ^' . , 

Jo s Jo s(s + l)™+p 

where we made the change of variable w = {s + l)u to obtain V^{x{s + l)p). Now 
we use Proposition ^ (iv) with a = (s + 1)^ > 1 to obtain 



1 r^ u"^\nue-- ^ r.,p. ./e'^ 1 

nm + l)Jo (xP + uf-'p - Jo '"^'^^U s{s + ir^^^ 

= V^{x)ij{m + 1) 



ds 



13 



when p > 1. For p < 1, Proposition ^ (iv) gives the inequahty in the opposite 
direction. Hence inserting the result in (|2lD yields 



dm 



Vmix) 



< if p > 1 
> iip> 1. 



To prove (ii) it is slightly more convenient to consider the logarithmic derivative 
ln[mV^(x)] and show that it is positive for p > 1 and negative for p < 1. 



dm 



Proceeding as above, we find for p > 1 

-— In mVZ(x) = h ^^ 

dm m 



< 



1 
m 

1 
m, 



-ds 



Vl{x) 
1 



- — ds — 

S .10 



m\ 

oo 



s{s + l) 



m+- 



ds — ipijn + 1) 



s[s 



1)' 



-ds — ipijn + 1) 



h ipim) — ipim + 1) 

1 roo 

m, 



\m+l 







'0 (s + iy 

where we have used (p2D and the following inequality with a 



[s + l) 



i/p 



V^iax) 



< 
> 



aVP(x) for 



p> 1 
p < I 



(23) 



for all a > 1. This is easily verified and implies that the inequality proved above 
for ^ In [mV^(a;)] is reversed when p < 1. QED 

The following result gives the asymptotic behavior of VJ^{x) for large x. 

Proposition 7 For p > 1, V^{x) has the asymptotic expansion 



(p-l)(m+l) (2p2-3p + l)(m2 + 3m + 2) 



r2p-l 



+ 



Proof: This follows from (|^) since 



2p2 X^ 



3p-l 



+ 



4p-l 



X' 



v:^{x) 



1 



r(m + 1) xP-i 

I I 

T{m + 1) xP-^ Jo 



° (l + Jr) ^ 



JEldu 



M™e-" 



\ {p-l)u ^ {p-l){2p-l)u^ 



xP- 



(p-l) r(m + 2) 
p r(m + 1) xP 



+ 



p x^ 2p^ x'^P 

(2p2-3p+i) r(m + 3) 
(2p2) r(m + 1) x^P ^ 



+ 



(iit 



2;3p 
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3 Recursion Relations and their Consequences 

3.1 Recursion relations for V^ 

Although the case p = 2 is of primary interest in apphcations, we continue to 
study general p is this section, as the proofs for general p are identical to those 
for p = 2. In these recursions, our convention that V^i{x) = x^~^ plays an 
important role. 



Proposition 8 For all m & H, m > 1, for all x > 0, 



V^ix) = - 
m 



p 



(24) 



Proof: For m = 1, one gets that 
,1 



Vlix) = pe' 



X- 



oo 

p\ I „-tp 



p 



e dt + X 



p 



xP)Vo^{x)+xP-^V^^{x). 



For m > 1, using (H) and integration by parts, we find 



VZix) 



pe 



r(m + l) 



mT(m 



if - x^)'^-^{f - x^)e-''' dt 



.-x'^ / e 

' X 



+ _ / e-*'((tP - x^)"^-^ + (m - l)ptP(tP - xP)"'"H^)dt 

PJx 

P Jx 



mT{m) 



+ (m - l)xP / {tP - xPr~^e-''dt 



{m-l + -- x^)Vl,,{x) + +xWl_,{x) 
V 



1 
m 

Repeated application of (^4]) gives a useful corollary. For m G R, let \m\ 
denote the "floor" of m, i.e., the largest natural number less than or equal to m. 

Corollary 9 Let m E'R., m>l and let n &'N such that n < \m\ . Then 

1 



efx) 



pm I 



{l-pxnVl_,{x) + VZ_,{x) + ... (25) 

+ . . . [p{m -n) + l]V^_^{x) + pxPV^^^_^{x) 
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In particular, if m is a positive integer, then 

m-2 

;i -pxnV^_,{x) + J2 V^{x)+px^V^_, 



pm 



^x) 

k=0 



(26) 



The expression (pB]) is well-defined for a; = 0. Putting a; = and using 
Proposition ^ (i), we obtain the (presumably well-known) identity 

r(m+i) 1 !^ir(A; + i) 

Pi = _ V ^ ^ r27) 

r(m + l) pm^r(A; + l) ^ ^ 

3.2 Averaged potentials 

These recursion relations are quite useful for studying the average of the first N 
of the Vm- For N a positive integer, we extend (^ to 

1 ^"' 
V^-'^ix) = -Yl ^™(^) (28) 

Tn=0 

Note that for p = 1, V^^^{x) = 1, for all x > 0. 

The next result follows immediately from (P^). 
Corollary 10 V^f{x) = pV^ix) - ^ [^(x) - V^_,ix)] . 
For the important case p = 2, this reduces to 

K^(x) = 2Viv(x) - — [V^,{x) - VW-i(x)] (29) 

The function Vo(|x|) is convex on (0, oo) but has a cusp at x = 0. How- 
ever, as discussed in property (f), for higher m both the convexity and cusp are 
lost. Thus, for higher m, the Vm are somewhat smoother than one might want 
for one- dimensional approximations to the Coulomb potential. The next result, 
although straightforward, is important because it implies that the averaged po- 
tentials Vj^{x) retain the cusp and convexity properties of Vq near the origin. 

Proposition 11 The function Vj^{x) is convex for all x > and 

hm ^Kv (a^) = -If- 
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Proof: Using (^) and (|T5D one finds 



— iVV;^(x) = Y, 2x[Kn(x) - Vm-i] = 2x[VW-i - V-i] = 2xVW-i - 2. 

m=0 



Therefore, to show that V^^ is convex, we need to show that 



IS uuiivcA, we 



xVat-i 



r(A^)Vo (l + :§)^ 



is increasing. This holds as, for x > 0, the function /i„(x) = [1 + u/x"^] ^^^ is 
increasing. 

Similarly, one can show that for p > 1, Vf^^ is convex on (0, oo) and 

.'}^+dx ^^ ^ ' N- 

For p < 1, the derviative becomes infinite at the origin; however, concavity of 
VP;^ on (0, oo) still holds. 

3.3 Polynomials defined by recursion 



We now observe that by repeatedly using ( [2^ ) to eliminate the V^ with the largest 
value of m from (p5|) allows us to write V^ in terms of the two "lowest" functions 
(e.g., Vq and V^i in the case of integer m) and that the coefficients in such 
expressions define two classes of polynomials related to confiuent hypergeometric 
functions. We discuss the properties of these polynomials in some detail. First, 
we make the statement above explicit. 

Corollary 12 Form > 1 there are polynomials P^iv) ^^^ Qmiv) of degree \rn\ 
such that 

Vlix) = P^{xnVZ_^^^{x)+xPQl_,{xnVl_^^^_,{x). (30) 

In the case of integer m (|30|) becomes 

Vlix) = P^{xnV,'{x)+xPQl^,ixnVl,{x) 

= P^ixnViix) + xQl,^,{xn (31) 

where the second expression follows from our convention V^i{x) = x^'^ . We 
define P^ (y) = 1 for m G [0, 1) and Ql,{y) = for m 6 [-1, 0). Then (0) holds 
trivially for ra G [0, 1). 
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Proof: The desired polynomials are defined recursively. First let 

pp (l/) = —fm-l + -~y] for m G [1, 2), and (32) 



Q'miy) 



m + 



— for m G [0, 1). 



(33) 



Then ( PP]) holds because it is equivalent to (^) for m G [1,2). For m > 2 define 

my) by 



m 



m-l + --y]P^_,{y)+yP}:^_,{y) 



(34) 



and for m > 1 define Q^{y) by 

1 



Q'miy) 



m + 1 



m + --y] (5™_i(z/) + y Ql,-2{y) 



(35) 



It is straightforward to use induction to check that (|24D yields (^). QED 



We now restrict ourselves to m G N and study these polynomials in more 
detail. The first few polynomials are given in the following Table. 



m 



pp 



Ql 








1 




1 

2 


1 

2 


-p-y 


1 

3 






The following useful results, which hold for m > 1, are easily checked by 
induction. B{x,y) = y(x+) ^^ ^^^ Beta function. 



Q^(o) 



r(m + i) 



r(m + l)r(i) mB{m,l) 
Tim + l + l) 1 



r(m + 2)r(i) (m+l)5(m + l,i) 
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P^iy) = - 

m 



m~l 



and 



Q'miy) 



1 



m + 1 



=0 



j=0 



(36) 
(37) 



We now obtain two expressions for ^Kn(^)- First, observe that using ( |3TD in 
^) yields 



Differentiating (|31| ) yields after some simplifications 
d 



dx 



V^{x) = px^-' [{P^r (x^K^ix) + P^ix^Kix) - V^,{x)] 



+ ^' {Qi-i)' (xnvi.ix) + iQi.,{xnv^,{x) 

where (-P^) iv) denotes ^P^iy)- Equating these expressions yields 

-[{pir{xn+p:L-i{xn]Vo'{x) = 

'x^ {Qi_,)' {xn - p^(x^) + (i - xnQi_,{xn + x^Qi_,{xn] n(^ 

This provides motivation for the following 
Lemma 13 For m G N, ??i > 1, 

d 



{3t 



dy 



P'.^{y) 



-Pra-iiy)^ and 



d 



y-Ql,_M = P^(y)-(m+l)g^(y)+mQ^_,(y) 



(39) 
(40) 



Proof: We first prove ( |39D by induction. It can be verified for m = 1, 2 using the 
Table above. Then using (|3^), we find 

= [m-i + l-y) ^yPl^Av) - Pl-Av) + y^/'ra-2{y) + PL^2{y) 

= -Pl-iiy) -{m-2 + l-y) P^_,iy) - y P^^_M 
= -mP^_-^{y). 
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This implies that the coefficient of VJf in (0) is identically zero. Therefore the 
coefficient of V^i must also be identically zero. Substituting y = x^ gives 

y^Q'm-M = P::r{y) + (y - 1) Ql-M - yQ'm-2{y) 



pp^iy) - (m + l)Qay) + mQl_,{y) 



where we used (|5B|). QED 



Note that since the left side of pO|) is a polynomial of degree m — 1, this 



implies the coefficients of the y"^ terms in P^ and (m + l)Q^{y) are identical. In 
fact, one can use (|34D and (|35D to see that the leading terms of P^ is {—l)"^y"^/m\ 
and that of Q^, is (-l)™?/'"/(m + 1)! 

A set of polynomials {p„(x)}belongs to the class known as Appell polynomials 
if they satisfy j^Pn{x) = Pn~i{x). Therefore, (|39D implies that for each fixed 
p, the set {—)"^Pm{y) forms a family of Appell polynomials. 

One can use (^) in ( |5^ to replace -P^-i ^^^1 P^-2 by derivatives of P^ and 
obtain a second order differential equation satisfied by P^. This allows us to 
obtain a relationship between the polynomials P^ and confluent hypergeometric 
functions, which we denote iFi(a,7,?/). 

Theorem 14 Form G N, m > 1, P^{y) satisfies the differential equation 

y4>"{y) - (m-lH y] (p'iy) - m4){y) = 0. (41) 

Standard techniques show that (0) has a polynomial solution of the form 
0(l/) = YJk=Q^ky'' with hk = - fcffr_\,) ^fc-i, A; > 1, &o 7^ arbitrary, and a second 

solution of the form 0(y) = ^'^^QCky'^^'^^^^P with Ck = - ^^^^^^^illf k-u k > 1, 
Co 7^ arbitrary. Since -Pm(O) = ' — T^ we conclude that bk = -—^ j- 

m B{m, — ) kl {m—k)B{m—k, — ) 

and 

^ (-l)'=r(m+i-A;) 

PP (y) = V ^ - ^ ^ - V'' 

^r(A: + l)r(m+l-A:)r(l)^- 

The restriction that 1/p be non-integer in the next result is neither serious, 
nor unexpected, in view of Proposition ^ 

Corollary 15 Let p ^ - for n G N. Then 

my) = r,] 1, e-y iFi (l - i, 1 - i - m, y) (42) 

m B{m, -j \ p p y 
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Proof: Write (f){y) = e ^<f){y)- Then it follows from ( ^I]) that (p satisfies 

yl"{y) -(m-l + ^ + y) $'{y) " (l - ^) ^v) = 0- (43) 

which has the form of the differential equation satisfied by the confiuent hyperge- 
ometric function. Comparing the behavior of P^{y) near y = with that of the 
well-known solutions to (|43|) suffices to complete the proof. QED 



It is well-known that for real a and 7, iFi(a,7,|/) has at most finitely 
many zeros on the real line. Hence, the same holds for P^. In fact, we can show 
that P^ has no zeros when m is even and exactly one when m is odd. 

To show this, it is convenient to introduce the new variable z = - — y and 

write PJ^{y) = -P^ i- — y]. The first few of these polynomials are Pq{z) = 1, 



Pf{z) = z, and P^{z) 



\ -" + 5 



Lemma 16 Form G N, ?7i > 2, the polynomials P^i^) satisfy 



(i) P^iz) = - 
m 



{m-l + z)P^_,{z)+{--z]P^_,y 



fii) Pliz) 



R tM^) 



1 
m 



m-2 



iX]Pf(z) + zP^,(z) 
.^ i=o 



, and 



Pm-li 



Proof: Follows immediately from substitution in 



and 



Corollary 17 For m ^ N , m > , all coefficients in the polynomials P^{z) are 
positive. 

Proof: This follows immediately from the explicit expressions above for P^ when 
m = 0, 1 and part (ii) of Lemma [T6| . 

Proposition 18 Let m G N, ?7i > 1 

If m is even, P^{z) > 0. 
If m is odd, P^(-2) has exactly one root Zm- 
Moreover, the roots form a strictly decreasing sequence with —m + 1 < -2m < 0. 
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Proof: First note that Corollary |T^ implies that -P^(^) — ^^ all ;z > and 
lim,^f^ P-iLiz) = oo. Now we claim that 



771 V 



,. rin / \ I -OO if m IS odd 
lim i^£(2) = < .p 

oo II m IS even 

This can easily be verified by induction using part(i) of Lemma |16| (i) above. 

For m odd, P^iz) has at least one root Zm- We now prove by induction that 
Zm is the only root if m is odd and that P^{z) > for all 2; if m is even. The 
induction hypothesis is easily seen to hold for m = 0, 1. Suppose it is true up to 
m—1 and consider P^, with m odd. Then, by Lemma p!6|(iii), ■j^P^{z) = P^_i{z). 
Since m — 1 is even, by the induction hypothesis P^_i{z) > 0. Thus P^i^) is 
increasing for all 2; G R, which implies that for m odd P^i^) has only one root. 
Since P^{z) > when 2; > that one root must satisfy z^ < 0. 

If m is even, then ^_P^(z) = P^_]^(2;) which, by the induction hypothesis, 
has exactly one root z^-i < 0. Therefore P^{z) has a local extremum at Zm-i- 
Since jj^P^{z) = Pm-2i^) ^^^ m — 2 is even, P^_2iz) > by the induction 
hypothesis and Zm-i is a local minimum for P^{z). By Lemma |16| (i) we have 



"^-Pm(^m-l) - \p~ ^m-lj Pm-2i^m-l] 



since P^_i{zm-i) = 0. But by the induction hypothesis 2;m-i < andP^_2(-2m-i) > 
so that Pf^{zm-i) > as required. 

It remains to be shown that the roots are decreasing and bounded below by 
— (m — 1). Both can be easily checked for m = 1, 3 and then proved by induction 
using Lemma |l^ (i). We now let m be odd. Since P^iz) is increasing to show 
that Zm > —m + 1, it suffices to show that P^{—m + 1) < 0. For z = —m + 1 the 
recursion relation reduces to 

mP^^{-m + l) = (l + m-l) P^_,{-m + 1) 

which is negative by the induction assumption hypothesis that -2^,-2 > ~fn + 3. 
To show that Zm < Zm-2 it suffices to show that P^{zm-2) > 0. But 

mP^{Zm-2) = (m - 1 + Zm-2)Prn-li^m-2) > 

since Pm-2i^rn-2) = 0, Zm-2 > —rn + 3 > —m + 1, and Pm-i{z) is positive. QED 
We now restate the results above in terms of the behavior of the original 
polynomials Pmiv)- 
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Corollary 19 Let m G N, -m > 1. Then 

(i) If m is even, then P^iy) ^ for all y G R. 
(ii) If m is odd, then P^{y) has exactly one root ym^O- 
(iii) For all m > 1, 

Urn Pl{y) = oo 

— oo if m is odd 



y~*oo I OO if m IS even 

Although we were able to obtain an explicit expression for the polynomials 
P^(x) relating them to confluent hypergeo metric functions and analyze their 
behavior in some detail, we do not much information about Q^{x). This is, 
at least in part, because (|40|) mixes Q^{x) and P^i^) and does not lead directly 
to a differential equation for Q'^{x). It would be interesting to know more about 
the polynomials Q^{x). 

4 Inequalities and Convexity 

4.1 Inequalities for Vo{x) 

We flrst illustrate our strategy by proving a special class of inequalities for Vq. The 
convexity of l/Vo(a;) follows directly from the optimal upper bound in this class 
as given in Theorem ^ below. Although, as discussed at the end of this section, 
these inequalities generalize to Vm the resulting upper bound is not sufficient 
to establish the convexity of l/Kn- For this we need a bound on the ration 
Vm{x)/Vm-i{x). Nevertheless, these simple inequalities for Vq, which can also be 
interpreted as ratio bounds, are of some interest in their own right in a variety 
of applications. Because the geometric strategy is also used in our more complex 
proofs of ratio bounds, we think there is some merit in presenting it flrst here. 
We now deflne 

9kix) = -—--—y===. (44) 

[k — l)x + va^ + k 

Theorem 20 For x > 

gA^) < K)(a;) < 9iix) (45) 
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and these inequalities are optimal for functions of the form (^^ with equality only 
at (7.(0) = 1^0 (0) = v^. 

Proof: It is easy to see that the family of functions gk{x) is increasing in k and 
that < (y'fc(x) < 1/x. In order to prove that the upper bound is optimal, we first 
observe that g'j.{x) = —klgkixyf'lx + {k — l)\/x'^ + k ]/\/x'^ + k and xgk{x) — 1 = 
—kgk{x)/[x + y/x"^ + k ]. Then one can verify that 

g'f,{x) > 2[xgkix) - 1] 



k {k-l)Vx^+k + x 2k 



\lx'^ + k {k — l)x + \/x^ + k X + \/x'^ + k 
^^ {k - 2)x^ + k{k - 3) < (A; - 2)xVx^ + k 
^^ x'^{k-2){k-A) + k{k-3f <0 (46) 

when k > 3. We now restrict attention to 3 < A; < 4 and let hk{x) = gk{x) — VQ{x). 
For fc = 4, the expression (^) implies g'^^x) < 2[xgk{x) — 1] so that h'^{x) < 

2xhi{x) for all X > 0; whereas for k < A this holds only for x < a^ = \/ (fc-2V'4-fc'i • 
Since both Vq{x) and gk{x) are positive and bounded above by 1/x, their difference 
also satisfies \hk{x)\ < 1 / x ^>- Q . 

For fc = 4, if hi{x) < for some a; > 0, then h'^{x) < 2xhi{x) is negative and 
thus hi is negative and strictly decreasing from a certain x on, which contradicts 
lima-^oo ^^4(3;) = 0. Thus h^i^x) > so that g^^x) > Vq{x), for all x. Now suppose 
that for some k < 4, g^ is an upper bound, i.e. hk{x) > for all x > 0. In 
particular, hk{x) > for all x > a^- For A; < 4, we find however that /^^(x) > 
2xhk{x) holds for x > au- Thus we get hk{x) > and strictly increasing for 
all X > Qk which contradicts lima;^oo hk{x) = 0. Thus the upper bound can not 
hold when x > at and k < 4. The lower bound also fails for k > it since then 
hk{0) = gkiO) - l^o(O) = Vk-y^>0. 

To establish the improved lower bound g.,r < Vo{x) we note that the argument 
above implies that hk{x) is negative for x > ak and 3 < fc < vr. However for 
k < 71 we have hk{0) < so that hk{x) is also negative for very small x. If hk{x) 
is ever non-negative, we can let b denote the first place hk{x) touches or crosses 
the X-axis, i.e., hk{b) = and hk{x) < for x < b. Then hk must be increasing 
on some interval of the form {xo,b). However, by the remarks above, hk{b) = 
implies 6 < a^ so that h'f^{x) < 2xhk{x) < on (xo,6). Since this contradicts hk 
increasing on (xq, 6), we must have hk{x) < for all x > if A; < vr. 

Thus we have proved the lower bound gk{x) < Vo(x) on [0, 00) for k < n. 
Since gk is continuous and increasing in k, it follows that g^ix) < Vo(x). To show 
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that this inequahty is strict except at a; = 0, note that the right derivative of hk 
at satisfies h',^{0) = 2 — k so that h'^{0) < and hT^{x) is negative at least on 
some small interval (0, xi). Then we can repeat the argument above to show that 
K{x) < if X > 0. QED 

As discussed in [^, |l^, [18| the upper bound implies the convexity of l/Vo{x) on 



(0, oo); in fact, it is not hard to use the fact that (|l5|) reduces to ^Vo(x) = 2(xVo — 
1) to see that the upper bound is equivalent to convexity. It was established 
independently by Wirth [Q and by Szarek and Werner fl^. (The latter actually 
proved slightly more by using (|1|) to define an asymmetric extension of Vo{x) to 
negative x. They showed in [|l^ that this extension is convex for x > 
Both bounds in (^5]) are sharper than the inequalities of Komatsu 



V2-' 



The weaker lower bound g^ix) < Vq{x) was used in ^ to show that the function 
[l/Vo(a;) — x] /Vo{x) is decreasing for x > 0. The lower bound g-rr{x) < Vo(x) 
was established earlier by Boyd as the optimal bound in a different class of 
inequalities. There is an extensive literature (see. e.g. fl^) on bounds for Vo{x); 
however, the class of inequalities obtained using functions of the form gk{x) does 
not seem to have been considered before so that the optimality of bounds of this 
type for k = tt and k = 4 seems new. 



Mascioni [|T2| generalized the upper bound to p > 2 for which he showed 
V,'{x) < ^^ 



3pxP-^ + a/p^x^p-s + 8p{p - l)xP-2 
and also showed that this implies convexity of 1/Vq{x) for p> 2. 



In view of Property (a) of Section |1.2| , it would seem natural to try to generalize 
([20| ) using functions of the form 



9Ti^) 



k 



(k — \)x + a/x^ + m + /c 



(47) 



Note that the functions g^ are increasing in k and that Xwa^k^oag^ix) 
Therefore Property (a) implies that 

g'^ix) < V^x) < hm g^{x). 

K— >00 

As gj^ is continuous in k, there must exist im and jm such that 

gZ{x)<VUx)<gTJx). 



(4J 
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However, we have not obtained explicit expressions for im and jm- One might 
expect that the optimal lower bound occurs when im is chosen to satisfy g^ (0) = 
Vm{0). However, numerical evidence shows that this is false; in fact, this choice 
for im does not even yield an inequality. 

4.2 Ratio Bounds 

One of our main goals is to show that the function ,, I , is convex for integer 



m > 1. The key to this is the realization that (|45| ) can also be rewritten to give 
bounds on the ratio Vo(a;)/Vli(x) = xVo{x). We now let 

[k — l)y — m + a/ (y + my + ky 
and note that xgk{x) = Gl{x^) so that ( ^5[ ) is equivalent to 

GUx') < xV,{x) = ^ < Gl{x% 

For integer m > 0, convexity of y -*".■, can be shown to be equivalent to 

<'m—l \X ) 



In addition to this upper bound, we can show the following 
Theorem 21 Let m G N, -m > 0. Then the inequalities 

G^'-\x^) < R^{x) < Gf{x^) (50) 

hold and are optimal in k for the class of functions of the form G'^{x^). 

The upper bound is optimal in k for all m. The lower bound is optimal in 
the sense that 8 is the largest integer for which the lower bound in ([50|) holds for 
all m. However, as we discuss at the end of Section |5.3| , for fixed m one can find 
k{m) such that G™(~^(a;^) < Rm{x) holds with k{m) > 8. 

Since G^{y) is increasing in both m and k, its behavior at zero and infinity al- 
lows us to also draw some conclusions about the optimality in m of (^Up. -Rm(O) = 
1-2^ and G^(0) = l-j^ for all k. Therefore, G'^(O) < /2^(0) < G^,(0) implies 
z/ < 771 — I and /i > m — I for all A;, k' . Thus, if we insist that m be integer, there 
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is no choice of k which allows m — 1 to be replaced by m in the lower bound when 

m > 0, or m hj m — 1 in the upper bound. This argument does not, however, 

1 1 

rule out the possibility of bounds of the form Gf, ^ (x^) < Rm{x) < Gy ^ (x^). 

To examine the behavior at infinity, note that 



1 4m + fc + 2,_/l 

1 4m + 6 , ^ / 1 
2^ + ^^ 



G^{y) = 1--+ I, +^71)^ a^d 



Rrr^i^) = l-- + ^7^ + 3 



where the asymptotic expansion for R^ follows from Proposition |^. It then follows 
that Rm{\/y) < Gk'(y) ™plies ;U > m + 1 — |. Thus m is optimal for the upper 
bound if A; < 4 and any attempt to decrease m would require an increase in k. 
Furthermore, fi = m implies A; > 4 so that the upper bound in (|50| ) is optimal in 

k. 



We postpone the proof of Theorem 21, which requires a lengthy computation 



even for the case p = 2, to the next section. Our proof uses induction on m. 
Therefore, we are able to establish (|50|) and the theorems in the next section 
only for m a positive integer. We believe that they are also true for non-integer 
m. However, a proof would require either a different method or independent 
verification of the upper bound for an initial range, such as — 1 < ?7i < 0. 

The ratio Rm{x) is of interest in its own right, and our results are sufficient to 
establish that it is increasing in x on (0, cxd). This is proved in the next section 



after Theorem ^, which uses a similar argument 

I 



Theorem 22 For m G N, the ratio Rm+i{x) = y'^(V is increasing in x. 



4.3 Convexity of 1/V^ 



We now prove some important consequences of Theorem ^ The first is 
Theorem 23 For all m G N, the function 1/Vm{x) is convex on [0, oo). 
Proof: We need to show that 

1 Y _2[Vmix)f-Vmix)iVmix))" 



Vm(x) J Vmix 



\3 



>0. (51) 



It follows from the differential equation (^) and the recursion relation ( ^^ that 
Vmix)" = 2 [Kn(x)(l + 2m + 2x^) - 2Vm-iix){x'^ + m)] 
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so that ( |5lD holds if and only if 

[V^{x)f{l + 2m - 2x2) + 2V;„_i(x)Kn(x)(3x2 - m) - Ax^[V^^^{x)f < 0. 

After division by [Kn-i(a;)]^ this can be rewritten as P[Rm{x)] < where 

P{z) = z^{l + 2m - 2x^) + 2z{3x'^ - m) - Ax"^. 

c 



-AC^ 



Writing the roots of P{z) = Az + 2Bz + C in the non-standard form ^ 
we find that C^lx"^) is either the smaller of two positive roots (when x"^ > m + ^) 
or the only positive root (when x^ < m + ^). Since P(0) < 0, in both cases, we 
can conclude that 

z < C^ix^) implies P{z) < 0. 

Therefore, it follows from the upper bound in Theorem ^ that P[Rm{x)] < 0, 
hence (|l]) holds. QED 

Proof of Theorem [22| : Using (|I3|) one finds that 



d 
dx 



Rm+i{x) =2x 



Rm+l{x) 
Rmix) 



- 1 



After rewriting this in terms of Vm and then using the recursion relation (^) with 
p = 2 to eliminate Vm+i, one finds that R'^^^{x) > if and only if 

2(m + l)[R^{x)Y - (2m + 1 - 2x'^)R^{x) - 2x^ < 0. 

The polynomial P{z) = 2{m + l)z'^ — {2m + 1 — 2x'^)z — 2x^ has one positive and 
one negative root, and K^^^ix) > if and only if Rmix) lies between these two 
roots. Since 1 > Rmix) > 0, it follows that Rm+\ix) is increasing if and only if 
Rmix) is less than the larger root, i.e.. 



Rmix) < 



Ax^ 



4(x2 + m)2 + 1 + 4m + 12x2 + 2x2 - 2m - 1 



c 



'B+VB'^-AC' 



where, we have again written the root in the non-standard form 

Then using the upper bound of Theorem ^, we see that it suffices to show that 

4x2 

Rmix) < 



< 



^ix^ + m)2 + 4x2 + 3x2 _ j^ 

4x2 



^4(x2 + m)2 + 1 + 4m + 12x2 + 2x2 - 2m 
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or, equivalently that 

y (x^ + my + 4x^ + 3x^ — m 



> a/4(x2 + m)2 + 1 + 4m + 12x2 + 2x^ - 2m - 1 
which is easily checked. 

5 Proof of Ratio Bounds 

The proofs in this section, ahhough elementary, are quite long and tedious. The 
details were checked using Mathematica. 

5.1 Differential inequality 

In order to prove Theorem |21], it suffices to establish the following 

Lemma 24 Let G™ be given by 1^4^). Then 
(V,For,n>l. |-G™(.^) < 2.(|SgL _ i). 

(ii) For m > 4, -G-(x^) > 2x(^^|^ - 1), 

but the inequality (ii) does not hold for m < 4. 

Proof: The proof is based on the elementary principle that if a function on the 
half-line is zero at the origin and increasing, then it is non-negative. Unfortu- 
nately, the actual verification is rather tedious and requires the repeated use of 
this principle. For simplicity, we put x^ = y and assume |/ > 0. Then (i) is 
equivalent to 

E„(,).(|g|-l)-l-Gr(,)>0. (52) 

Let Bm = {rri^ + y'^ + Ay + 2my) 3 . Then 

GTiy) ^^ 



Bm + 3y -m 
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and 

B^ [Am + {Bm + 3y - m)(E^_i - Bm + 1)] - (4m^ + 8y + 4my) 

Thus Em{y) > if and only if 

Bm ( 3m + Bm-i [Bm + 3?/ — m] j + m^ + 2my 
> 3m^ + 4y + lly^ + m^y + 5my^ + 3y^ + 5^ ('(y + nif + y") . 

We put 

/ \ O o *? 

s = s{y,m) = 3m + Ay + lly + m y + 5my + 3y , 
t = t{y,m) = 2my + 171"^, and 
h = h{y,m) = {y + m)'^ + y — 3m. 

Then E^iy) > if and only if 

BmBra-i{Bm + 3y - m) > B^h + s - t. (53) 

Notice that both sides of (^) are positive. For the left side this follows immedi- 
ately from Bm > m. For the right, note that Bmh{0) + s(0) — t(0) = and 

— [Bmh{y) + siy)tiy)] 

= —— [6y + Vlmy + 9m^y + 12y'^ + Qmy^ + 3y^ 

+ Bm{22y + lOmy + 9?/^) + 45^ - 2mBra + 3m^ - 6m + m^Bm] ■ 

Now observe that 

ABm — 2m.Bm + 3m^ — 6m + m.'^Bm 

= 3m{m^ - 2) + ^^(m^ - 2m + 4) > 3m(m2 - 1) > 0. 

since m > 1 and -B^ > m. Hence -Bm^ + s — t is increasing in y and the right side 
of ( |53| ) is also positive. Therefore we can square both sides of (|3]) to conclude 
that it is equivalent to 

Fiy)=Bmf,iy)-f2iy)>0, (54) 
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where 

fi{y) = {rn + y){y^ + my'^ + 3y^ — m^y — 3'my + 2y — m^ + 2m^) 

and 

f2{y) = y^ + "imy^ + by^ + 2m^y^ + bmy'^ + 2y'^ 

— 2m^y'^ — Sm^y^ — 3m^y — m^y + Qm^y — rrt' + 2rn^. 

Note that F(0) = 0. Therefore, to prove ( |5lD it is enough to show that ^F(y) 
Bn.f[{y) + ^^^^^^r^^"^ - n{y) > 0, or equivalents 



D{y)^d,{y)~BrrMy)>0, (55) 



where 



d^{y) = Blf[iy) + My){2 + y + m) 

= 6r7i^ — m'^ — 3m^ + 12my + 4?7i^y — 13?7i^?/ — Tm^y + 20?/^ 
+ lAmy'^ + Tm^y^ + 2m^?/^ + 48?/^ + 49my^ + 18m^y^ + SOy"^ 
+ 17my'^ + 5y^, and 
c^2(y) = f2iy) — 6"^^ — "^^ — 3m^ — 6?7i^|/ — Arn?y 

+ 6y^ + 15my^ + Qrn^y'^ + 20?/^ + I2my^ + 5?/"^. 

Note that -D(O) = 0. Therefore, to prove (|55D it is enough to show that j-D(y) 
d\{y) - B^d'^{y) - ■^^(^)(2+?/+"^) > 0, or equivalent^. 



G{y)=B„My)-92{y)>^. (56) 



where 



9i{y) = d[{y) 

= 12m + Am^ - 13m^ - 7m^ + 40?/ + 28m?/ + lAm^y + 4m^?/ 

+ 144?/^ + lA7my'^ + 54m^|/^ + 120?/^ + 68m2/^ + 25?/^, and 
92{y) = Bld'^{y)+d2iy)i2 + y + m) 

= 12m? + Am? — IStti^ — 7m? — ISm^y — l^m^y — 3m.'^y 
+ 602/2 + ISOmy^ + 183m^y^ + 58m^y^ + 298y^ 
+ 353my^ + 122m^y^ + ITOy^ + 93my^ + 25y^ 
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Note that G{0) = 0. Therefore, to prove ( |56D it is enough to show that j-G{y) = 
Brng'M + ^^^y^^y"'^ - g',{y) > 0, or equivalently 

H{y) = h{y)~B,nh2{y)>0, (57) 

where 

hi{y) = Blg[{y) + g^{y){2 + y + m) 

= 24m + eOm^ + Qm^ - 13m^ - 3m^ + 2A0y + 300my + AQOm^y 

+ 3A7m^y + llSm^y + 1520y^ + 22AQmy^ + IGGSm'^y'^ + A82m^y^ 
+ 21 I2y^ + 2233my^ + 738m^y^ + 930y^ + 497m7/^ + I25y^, and 

h2{y) = g'2{y) 

= —ISm^ — 13m^ — 3m^ + 120y + 360my + 366m^y + llGm'^y 
+ 894?/2 + 1059m|/2 + 366m^y'^ + 680y^ + 372my^ + 125y\ 



Note that H{0) = 12m{2 + 5?7i + 2'm'^) > 0. Therefore, to prove ( |57D it is enough 
to show that j^H{y) = h[{y) - B^h'^ly) - ^^(^)(g+?^+'") > o, or equivalently 

h{y)Bm-l2iy)>0, (58) 

where 

h{y) = Kiy) 

= 240 + 300m + 460m2 + 347m3 + 113m^ + 3040?/ + 4492m?/ 
+ 3326m'^y + 9QAm^y + Q33Qy^ + QQ99my^ + 2214m\^ + 
+ 3720?/^ + 1988my^ + 625y\ and 
kiy) = Blh',iy) + h2iy)i2 + y + m) 

= 84m2 + 316m^ + 347m^ + 113m^ + 720?/ + 2520m?/ + 5046m^y 
+ 3899m^y + 1077m^?/ + 9180?/^ + 15780m?/2 + 11727m2?/^ 
+ 3178m^y^ + 12202?/^ + 13145m?/^ + 4202m\^ + 4970?/^ 
+ 2613my^ + 625?/^ 

Note that h{y) > and hiy) > for all y > 0. Therefore (^) holds, if and only 
if L{y) = i?^(/i(?/))^ — (^2(y))^ > 0, which follows immediately from the fact that 
all the coefficients are positive in 

L{y) = 4(14400m2 + 36000m^ + 75936m^ + 97368m^ + 78972m'' + 37188m^ 
+ 8136m^ + 57600?/ + 172800m?/ + 717360m2?/ + 1373400m^?/ 
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+ 1732428m^i/ + 13603Um^y + 599454m> + llUAAm^y + ISAAOOOy^ 

+ 3838560m|/2 + 8A37260m'^y^ + 11062920m^y2 ^ 8495031m^y2 

+ 3499083m^y2 ^ sgsggGm'^y^ ^ 9342880^^ + 24217360my2 

+ 32546720m^y^ + 24561680m^y^ + 9950080m^y^ + 1694280m^y^ 

+ 17918380?/^ + 37038224my^ + 34271234m\^ + 15627870m^|/^ 

+ 2862630m^y^ + 15343236?/^ + 23700982my^ + 13930330m\^ 

+ 2982516m3?/^ + 64459637/*^ + 6618363m|/^ + 1887294m22/^ 

+ 1302640?/^ + 667056m|/^ + 101250?/^). 

To prove (ii) we proceed similarly, but now let Bm{y) = \/{y + raY + 8?/ and 






E„(l/) = %G^iy) - (^^ - 1), noting that 



GTiy) 



Bm + 7y -m 

We now need to show that Em{y) > for all m > 4. As the argument is similar 
to that above, we omit the details except to indicate the steps leading to the 
condition m > 4. Observe that Em{y) > if and only if 

Bm{{y + mf + y -7m) 
+7m^ — m^ + 24y — 2my + 5m^y + 55y^ + 13my'^ + 7y^ 

>BmB^m~i{7y-m + Bm). (59) 

Again, both sides of the inequality are positive. Hence we can square both sides 
of the inequality and, as above get that (^) is equivalent to 

F{y) = fi{y)-B^f2{y)>0, (60) 

with the appropriate f\ and /2. Again -F(O) = 0. Therefore, in order to prove 
(|60|), it is enough to show that j-F{y) > 0, or equivalently, after rewriting, 

D{y) = Bmdi{y)-d2{y)>0, (61) 

with the appropriate di and ^2- And again D{0) = 0. We repeat the procedure: 
To prove (|6ll) , it is enough to show that j-D{y) > 0, or equivalently. 



ei{y) - B„,e2{y) > 0, 

with the appropriate ei and 62- ei and 62 turn out to be both positive for y > 0. 
Therefore (|6T|) holds if 

L{y) = {e,{y)f~{B^e2{y)f>Q. 
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L(0) = and 

L'(0) = 192m^(m - 4)(1 + 2m)(480 + 64m + 90m^ + 33m^). 

Thus L'(0) > if and only if m > 4. For all m > 4, L"{y) > for all y>0. This 
finishes (ii). QED 



5.2 Proof of Theorem 

We will prove by induction that Rm{x) < G^i^x^) for m = 0, 1, 2, 3, . . . As ob- 
served earlier, this inequality holds for m = 0, since it is then equivalent to the 
upper bound in (|^). Let 



HUx) = G^(x2) 



Rm{x). 



Then the upper bound in Theorem ^ is equivalent to Hm{x) > 0. This can be 



verified using the strategy of Section |4.1| if the following conditions hold 
(i) HmiO) > 0. 



fii) lim Hry^ix) 



0. 



iii) H'^{x) < Fpos{x)Hm{x) for some strictly positive function Fpos{x) > 



Conditions (i) and (ii) hold. Indeed, 



HmiO) 



2m 



T('m)T(m + ^ 



l + 2m r(m+l)r(m 



1 



2m(l + 2m)' 



and 



lim Hjn{x) 



0, 



since lim Rm{x) = 1, and for all /c > 1, lim G^{x) = 1. 



We check now condition fiii). It follows from Lemma 

Rm(x) 



(i) and (E% that 



H'(x) < 2x 



GTix^] 



G 



2x 



GT-Hx^)Rm-i(x 



Rm-lix 

[G'^{x^)Rm-i{x)-G 



4 \X )rCrn\X 



< 



2x 



R'^-^ix] 

2x 



R 



rn—l I 



x) 



[G'^ix') - Rm{x 

Hmix) 
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where the inequahty follows from the induction hypothesis Rm-i{x) < G™~^(x^). 
Thus (iii) holds with -Fpos(a;) = 2x/Rm-i{x). 

For m > 4, the lower bound is proved similarly. One considers Hm{x) = 
Rm{x) — G™(a;^) instead and uses Lemma [2^(ii). The cases i?i, R2 and R^ have 
to be verified directly as Lemma ^ (ii) only covers the cases Rm for m> A. 

Using (p4D, _Ri > Gg is equivalent to showing that 



X 9x + Ux^ + (2x2 - l)V8 + x2 



Vo{x) 2(7x + VsT 



x^ 



As Vo is always positive, this inequality holds trivially for those x, for which the 
right hand side is negative or zero. Therefore we only need to prove the inequality 
on the interval [xq, cxd), xq — 0.2511, where 



9x + 14x^ + (2x2 - l)V8 + x2 > 0. 
Hence we need to show that for all x G [xq, cc), 

T/ r w 9 7x + V8 + x2 

Vn X < 2x 



9x + 14x3 + (2x2 - l)v/8Tx2 



6x2 _^ 

2x ^=^= (62) 

1 + 6x2 + i2a;4 - 2x^8 + x2 ' 



Put hi{x) = 23; ^^g^,^^^^^j4_2^^8+x^ - -^y Theorem |2g of Section |A|, inequality (|6 
is true for all x G [xq, 00), for which 

5(4 (x) < hi{x). 

This last inequality holds only on an interval [xo,Xi], Xi ~ 1.399. For all x > Xi, 
we show that 

hi(x) < — 

X 

and 

K < 2(a;/ii - 1). 

Then ( |62D follows as in Section ^A\ 

Next, we find that R2 > Gg is equivalent to Vo(x) > h2{x), where 

, , , ^ 3 + 9x2 + 14x^ + (2x2-3)(8x2 + (l+x2)2))l 

-3 - 7x2 + 32a;4 ^ 28x6 + (3 - 4x2 + 4x4) (8x2 + (1 + x2)2))l 
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and i?3 > Gg is equivalent to Vq{x) < h^{x) where 

2x 

N{x) 



2x r 
/i3(x) = — -- -30 - 23x2 + 32x^ + 282;^ + ^/Sx^ + (2 + x^y (15 - Sx^ + 4x^ 



with 

N{x) = 30 + 3x2 _ 42a;4 ^ 923,6 ^ ^g^s 

+ (8x2 ^^2 + x2)2))5(_i5 + 18x2 - 12x^ + 8x^). 

Again, we have to check these inequahties for Vq only for those x, for which the 
right hand sides are positive. We then proceed as for Ri and show that g-j^ > /i2 
up to a certain X2 and that /i2 < ^, h'2 > 2{xhi — 1) on [x2,cx3). Similarly, we 
show that g4 < h^ up to a certain X3 and that h^ < -, h'^ < 2{xhi — 1) on [X3, 00). 
QED 

Note that these arguments also show that on the interval [xi, cxd) the function 
hi is a better upper bound for Vq than g^] on [x2, 00) the function /12 is a better 
lower bound for Vq than g^,] and on [xs, 00) the function h^, is a better upper 
bound for Vq than g^^. In fact, h^ < hi < g^ for x > X3. 

5.3 Optimality of bounds 

We still need to consider optimality of the lower bound in upper bound in (^ in 
the parameter k. We continue the strategy above using similar notation so that 

now B^ = ^{y + mf + ky and E^iy) = [^^ - l] - %GT{y) with 

hii 

GTiy) 



Bm + {k- l)y -m' 

Then E^iy) < if and only if, 

2B^Bra-i{B^ + {k- l)y - m) < 2B^{y + {y + mf - (fc - l)m) + P (63) 

where 

P = — 2m2 + 2km^ — 2m^ — 2ky + k'^y — Amy — Qnn?y + 2km'^y 
- 22/2 - 2A;y2 + 2A;2y2 _ Q^y^ + Akmy^ - 2y^ + 2ky^. 

For m > 1 and k > 2 both sides of the inequality are positive. Therefore we can 
square both sides and get that E^iy) < if and only if 

Fiy) = fiiy)-B^My)>0, 
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where 

f^(y) = -IQkm^ + Sk^m"^ + 32m^ - 16km^ - 20k'^m'^y + Sk^Tii^y + Skm^y 

- APm^y + 160m^y - 9Qkm^y + Sk'^m^y - Ak^y^ + k^y^ 

+ UOkm'^y^ - SAk'^m'^y^ + Uk^m^y^ + 320m^y'^ - 22Akm^y^ 
+ 32Pm^y'^ + 20Py^ - 20k'^y^ + Ak'^y^ + lh2kmy^ - 12APmy^ 
+ 24k^my^ + 320m'^y^ - 256km^y^ + 48k^m^y^ + 56%^ - 52Py^ 
+ 12kY + IGOmy^ - lUkmy^ + 32k^my^ + 32y^ - 32ky^ + Sk'^y'^, 
fiiy) = 4:[-Akm^ + 2Pm^ + 8rn^ - Akm'^ -3Pmy + k^my + 2km^y 

- k^rr?y + 32rrC'y - 20km^y + 2k'^m^y + k'^y'^ + IQkmy'^ - 12k'^my'^ 
+ 2k^my'^ + A^m^y"^ - 3Qkm^y'^ + Qk'^m^y'^ + 10%^ - dk'^y^ + 2k^y^ 
+ 32my^ - 28kmy^ + Gk'^my^ + 8y^ - 8ky^ + 2k'^y^) . 

Then -F(O) = and in order that F > 0, we must have ^F(O) > 0. Computing 

j-F{y), we find that j-F{0) = 0. We apply the same procedure as in the proof of 
Lemma ^, compute the successive derivatives and evaluate them at 0. Evaluating 
the derivative at 0, in the fourth step of the procedure gives the value 

24fc^r?7,(l + 2m){km — 6m — k). 

Therefore, in order that (^) (which is the condition for the lower bound) holds 
for all m > 2, we have to have at least that k > linim^oo ^^ty ~ ^" Thus for 
m = 2, k > 12 will do, for m = 3, k > 9, for m = 4, k > 8 and so fourth. 
Therefore, as G™ is increasing in k, it seems a natural choice to pick /c = 12 or 
bigger for the lower bound. And indeed, one can check that G^ satisfies the lower 



bound condition of Lemma E3 for A; > 12 and m >2. However, it is not true that 



for k > 8, G^ ^ is a lower bound for Rm, for all m > 1. It is a lower bound for all 
m > m{k), from a certain m{k) on. Thus the induction in the proof of Theorem 



21] cannot start at m = or m = 1. For m < m(k), there exists Xm such that 
Rm - G^-^ > on [0, x^] and Rm - G^'^ < on (a;„, oo). 

5.4 Extensions to general p 

For p = 1, all the functions involved are identically equal to 1 and hence trivially 
convex. For large x, 1/V^{x) ~ x^~^ and x^~^ is concave for 1 < p < 2. Hence 
we can not expect convexity of l/V,^ on (0, cxd) for p in (1,2). It was shown 
in [0 that yp is not convex on R^ for < p < 1. Therefore, we study only 



generalizations to p > 2. Our method of proof yields verification of the convexity 
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of Tjjr for all m > 1 up to at least p = A. However this method breaks down for 
larger p. 



We generahze our previous notation to R^i^) 






and observe that 



d 

dx 



Rtix) 



px 



p-i 



Rl^{x) 



T3P 

1 Lyy 



m—l 



X) 



(64) 



For /i;>l,m>0, p>lwe generalize G™ to 



G 



m,p . 



X'^ 



kpx^ 



p[{k — l)xP — m] + ^Jp^{xP + Wb)^ + 2kp{j) — 1)xP 



The proofs of Theorems BT^, 22 and E3^ can be extended provided that the 



analogue of Lemma ^ holds. This is not the case for large p. However, although 
the generalization of the upper bound in Theorem |2T] is a necessary and sufficient 
condition for convexity of 1/V^(x), Lemma ^ is only a sufficient condition for 



Theorem 2L Indeed, we were able to establish the lower bound in Theorem 21 for 



m 



1, 2, 3 even though part (ii) of Lemma |2^ does not hold for m < A. Hence 



the fact that Lemma ^ breaks down for large p does not preclude convexity of 



VKn(^)- O^ th^ contrary, numerical evidence suggests that -^ is convex for all 
p>2. 



Lemma 25 For all A> p > 2, m > 1, 

d 



dx 



GT'^ix^) < px 



p-i 



anxn 



G 



'in—l,p I 



X'P) 



-1 



This is equivalent to 



EL 



px 



p-i 



G 



'">-^P(rj.P 



G 



m—l,p. 



XP] 



d 

dx 



G 



m,p / 



>0 



which allows us to make some remarks about the range of validity. Although 



Lemma |2g can probably be extended to some higher p, it does not hold for all 
p, m. On the contrary, for all m > 1 there exists p{m) and an interval (x^ , Xg ) 
such that EP^ < on that interval for all p > p{m). For example, numerical results 
show that for ?7i = 1, an interval on which Em < exists when p > 10; for m = 2 
when p > 14 and for m = 3 when p > 18. 



For simplicity, we only sketch the proof of Lemma ^ and give the final ex- 
pressions for p = 3. Similar expressions can be given for p = 4. We have checked 
the details using Mathematica, but omit the long formulas. As the expressions 
involved are monotone in p, this suffices for the entire interval 2 < p < 4. 
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Proof: Let Bm = \/v^{y + '^Y + 8p(p — I)?/- Then E"^ > is equivalent to 

Brr,Brn-liB,n + Spy - pUl) 

> pBm [p{y + mY + y{hp — 8) — Spm\ + p^ [Sm^p — w?p 

+ |/(— 8 + 8m + 8j9 — Qmp + m^p) + |/^(— 24 + 23p + 5mp) + 3py^} . 

Following the procedure used to prove Lemma |2^ (i), we eventually find that it 
would suffice to show that 

/i(y)5„-/2(y)>0, (65) 

where, for p = 3, 

l^{y) = 3(5120 + 2880m +13800m2 + 11034m^ + 3051m^ + 169600y 

+ 199632my + llQ820m^y + 26028m^?/ + 2982962/^ + 239706my^ 
+ 59778m\2 ^ 13392O1/3 + 53676my^ + 16875?/^), and 

/2(y) = 27(-2048m + 192m^ + 3448m^ + 3678m^ + 1017m^25600|/ 
+ 49920my + 76152m'^y + 45330m^y + 9693m^y + 200000y^ 
+ 250152my2 + 139266m2y2 + 28602m^y^ + 195880y^ 
+ 157254my^ + 37818m^y^ + 59640y^ + 23517my^ + 5625?/^). 

Both h{y) and hiy) in (^) are positive, and hence (|^) is equivalent to 

L = h{yYBl-hiyf>0. 

which can be verified for m, > 1. QED 
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